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Abstract 

Numerical simulations of coupled map lattices (CMLs) and other complex 
model systems show an enormous phenomenological variety that is difficult to 
classify and understand. It is therefore desirable to establish analytical tools 
for exploring fundamental features of CMLs, such as their stability properties. 
Since CMLs can be considered as graphs, we apply methods of spectral graph 
theory to analyze their stability at locally unstable fixed points for different 
updating rules, different coupling scenarios, and different types of neighbor- 
hoods. Numerical studies are found to be in excellent agreement with our 
theoretical results. 
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1 Introduction 



Coupled map lattices (CMLs) are arrays of cells whose state value is continuous, 
usually within the unit interval, over discrete space and time. Starting with Tur- 
ing's seminal work on morphogenesis pQ , they have been used to study the behavior 
of complex spatiotemporal systems for more than 50 years. More recently, Kaneko 
and collaborators have established many interesting results for CMLs [2] as gener- 
alizations of cellular automata, whose state values are discrete. 

One key motivation for modeling spatiotemporally extended systems with CMLs 
is to simplify the standard approach in terms of partial differential equations. And, 
of course, CMLs would not have become accessible without the rapid development 
of fast computers with large storage capacities. Within the last decades, CMLs have 
been applied to the study of diverse as social systems, ecosystems, neural 

networks, spin lattices, Josephson junctions, multimode lasers, hydrodynamical tur- 
bulence, and others (cf. the special journal issues CHAOS 2(3), 1992, and Physica 
D 103, 1997). 

A compact characterization of a CML with one time parameter is given by 

u(n + 1, x) = (1 - e)f(u(n, x)) H V" g(u(n, y)) , (1) 

Tlx 

where x represents the sites of the lattice, x = l,...,N ta t, here considered as the 
vertices of a graph, and n represents the time step of the iteration. The parameter e 
specifies the coupling between each cell and its neighborhood (and is often considered 
as constant over time and space). The sum over y ~ x is the sum over all n x 
neighbors y of vertex x. The function g characterizes the interaction of a vertex 
with its neighborhood and will be explained below. 

As in many other studies of CMLs, f(x) is the logistic map on the unit interval, 

f(x) = rx{l — x) , 

with < r < 4. For r > 1 the logistic map has a critical point at — - which is 
unstable for 3 < r < 4. The relevance of maps with quadratic maximum (such as 
the logistic map) for models of neurobiological networks was recently substantiated 
by novel results concerning a non-monotonic (rather than sigmoid) transfer function 
for individual neurons |3 . 

For e —>■ 0, there is no coupling at all; hence, local neighborhoods have no 
influence on the behavior of the CML. This situation represents the limiting case of 
N to t independently operating local objects at each lattice site. In the general case 
< e < 1, the independence of individual cells is lost and the lattice behavior is 
governed by both local and global influences, depending on the chosen neighborhood. 
CMLs with a maximal neighborhood, n x « N tot , are often denoted as globally 
coupled maps. Their behavior is determined by global properties alone (mean field 
approach) . 

The function 

g ( x ) = ax + {3f(x) + 1 f(f(x)) , (2) 
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with 

a + /? + 7 = l a,/3,7>0, 

allows us to treat the interaction between each vertex and its neighborhood in dif- 
ferent ways, depending on its time scale At. If the interaction can be regarded as 
simultaneous, At « 0, the situation can be approximated by a = 7 = and (3=1. 
Such a type of coupling, sometimes called "future coupling" |E], will be referred 
to as non-causal coupling [3] in the following, since the simultaneity of the inter- 
action between vertex and neighbors makes a clear distinction of cause and effect 
impossible. 

The situation of a finite interaction time At > can be properly modeled by 
j3 = 7 = and a = 1. In this way, past states in the neighborhood of a vertex are 
considered to act on the present state of the vertex with limited signal speed, so 
that the effect of an interaction is delayed with respect to its cause. Such a type of 
coupling will therefore by denoted as causal coupling in the following. Corresponding 
lattice behavior has recently been studied in [U El HI El IE] • 

A third, somewhat exotic possibility arises for a = (3 = and 7 = 1. This case 
reflects the idea to model the action of future states of a vertex neighborhood on a 
present vertex state. More precisely, this refers to locally extrapolated future states 
and is justified for small e since then u(n + 1, y) ~ f(u(n, y)). In this interpretation, 
the case of non-vanishing 7 is in contradiction with causality; thus we refer to such 
a situation as anti-causal coupling. 

The inclusion of anti-causal coupling can be interesting if one wants to study the 
consequences of a decomposition of a fundamental time-reversal invariant evolution 
of a system into forward and backward components. Using such a decomposition, 
one can investigate the influence of both components on the stability properties of 
CMLs. It has indeed been found that the stability of CMLs is supported by large 
a, and their stability is obstructed if a is small Based on these observations, 
one may speculate that stability acts as a selection criterion for a forward arrow of 
time HOj. 

Another time scale important for the physical interpretation of eq. (1) is the time 
interval Ar assumed for the updating mechanism, i.e. for the physical integration 
of signals from the neighborhood states with the vertex state considered. If signals 
between vertices are transmitted much slower than the time scale assumed for the 
updating mechanism, At <C At, the updating can be implemented (almost) instan- 
taneously, or synchronously. If this is not the case, At > At, updating must be 
implemented asynchronously. This entails the additional problem of determining a 
proper updating sequence, which can be random or depend on particular features 
of the situation considered. 

Most of the work on CMLs published in this respect (cf. JI]) was based on syn- 
chronous updating. For asynchronous updating as, for instance, studied by Lumer 
& Nicolis jT^j, it was found that the behavior of CMLs differs strongly from that 
of CMLs with synchronous updating. Additional results for asynchronous updating 
were reported by [T3J EH El El E] Asynchronous updating rules have been suggested 
as particularly relevant for neurobiological networks. 
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As a common feature of the (so far) few studies of asynchronous updating, it 
has been reported that it facilitates the synchronization and stabilization of CMLs 
decisively. In particular, Mehta & Sinha demonstrated that the dynamics at indi- 
vidual lattice cells is strongly synchronized by coupling among cells. Atmanspacher 
& Scheingraber jSJ |§| showed that unstable fixed points at individual vertices can be 
stabilized as a consequence of their coupling to neighboring unstable fixed points. 

Such a stabilization is of particular interest since it is independent of external 
control mechanisms. The global stabilization of unstable local behavior operates 
inherently, without external adjustment, once the coupling is strong enough. Such 
a possibility represents a powerful alternative to external control procedures in the 
style of "controlling chaos" [TKj . 

Most published results concerning the stabilization and synchronization of CMLs 
have been obtained by numerical studies for various parameters due to different 
coupling, different neighborhoods, and different updating. A number of theoreti- 
cal approaches have been developed to describe these numerical results analytically 
|16| IT7| ITHl \T§\ 120] ■ I n this contribution we present, for the first time, a compre- 
hensive stability analysis of CMLs at locally unstable fixed points, which allows 
us to understand corresponding numerical results theoretically. In particular, the 
sophisticated dependence of the stabilization behavior on the different parameters 
mentioned above jSJ |U] can be clarified. 

As recently proposed by Jost & Joy [21] and by Belykh et al. [22], the stability 
properties of CMLs can be compactly and conveniently analyzed in a graph theo- 
retical manner. We first determine, in section |2l the eigenvalues of the adjacency 
matrices related to various neighborhoods on regular lattices. In section |31 we derive 
a stability condition for synchronous update and constant solutions of the logistic 
map and apply it to special cases for causal and non-causal coupling numerically 
studied in jSJ E] . In addition, theoretical and numerical results are obtained for the 
situation of anti-causal coupling. In all cases, the agreement between theoretical and 
numerical results is excellent. In section 01 we present an intuitive argument for a 
stability analysis with asynchronous update. In section the discussion is extended 
to generalized neighborhoods. Section |U] summarizes and concludes the paper, and 
some perspectives are addressed. 

2 The Spectrum of a Graph 

A CML can be considered within the more general framework of the theory of graphs, 
whose properties can be studied by well-developed analytic methods (221 121] • A 
graph is defined by a set F of vertices, a set E of edges, and a map 8 ± : E — > F 
indicating whether a vertex is an initial point or a terminal point for an edge. A 
CML can be implemented on a graph with no self-loops, no multiple edges, and only 
undirected edges. Such graphs are sometimes referred to as simple graphs (cf. |25j). 
Many interesting properties of CMLs, in particular their stability properties, can 
then be investigated by two basic matrices characterizing the corresponding graph: 
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its adjacency matrix and its valence matrix. 

For an arbitrary simple graph (defined by the "neighbor-relation" ~ between its 
vertices), the adjacency matrix A(x,y) and the valence matrix V(x,y) are defined 
by: 

1 for y ~ x 
otherwise 



A(x,y) 



and 



V(x, y) = [l* with n x = J2 A(x, y) 



where n x is the number of neighbors of vertex x. The matrix A contains only Is 
and Os (no multiple edges), is non- reflexive (zero diagonal since there are no self- 
loops) and symmetric (undirected edges). If the valence matrix is proportional to 
the identity matrix, V — v 1, the graph is called regular. 

For the eigenvalues A of V~ X A we have |A| < 1 since V~ l A is a Markov matrix, 
i.e. {y~ x A)ij > and J2j(V = 1- The eigenvalue A max = +1 is assumed 
for the constant vector u(x). However, we will see that the lower bound for the 
eigenvalues of V~ 1 A can deviate remarkably from —1. From XxV~ l A = we can 
only deduce that the lower bound has to be negative (for infinite graphs it can be 
zero). 

In the following, we consider regular graphs, V — v 1. We determine the eigenval- 
ues for the adjacency matrix for different types of neighborhoods on a 2-dimensional 
N x N square lattice with periodic boundary conditions. In particular, we focus 
on the following types of neighborhoods: von Neumann 1st order (the four nearest 
neighbors horizontally and vertically), Moore 1st order (the eight nearest neighbors 
including the diagonal neighbors), von Neumann 2nd order (Moore 1st order plus 
the next to nearest neighbors horizontally and vertically), and Moore 2nd order (the 
24 lattice sites in a 5 x 5 square around the central point). 

The eigenvalues are calculated by introducing the shift matrix: 

/ 1 ••• \ 
1 ••• 



p = 



••• 1 
y 1 ••• J 

P is a unitary matrix with P N = 1 and P^ 1 = P + = P N ~ 1 . Therefore, the 
eigenvalues u k of P satisfy = 1 and are given by the N square roots of 1: 

(27riA; \ 
— j with k = 0,1, . . . , N — 1 . 

On a 2-dimensional N x N lattice we define the two shift operators for each direction 
by: 

Pi = P (8) 1 and P 2 = 1 ® P . 
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The adjacency matrices for the different neighborhoods commute with Pi and P2 
and can be completely expressed in terms of these operators: 

A v .N.l = Pl+Pf 1 + P 2 + P 2 - 1 

A MA = (i + p 1 + p r 1 )(i + p 2 + p- 1 )-i 

Ay. N. 2 = (1 + Pi + Pf 1 ) (1 + P 2 + P 2 V ) ~ 1 + Pi + Pf 2 + P| + P 2 ^ 2 

A M . 2 = (1 + P 1 + Pf 1 + P 2 + Pf 2 ) (1 + P 2 + P' 1 + Pi + P- 2 ) - 1 . 
Therefore, the eigenvalues of A are given by: 
v. N.l: 2 cos 0i + 2 cos 02 

2irki 

with (pi = — where fa = 0, 1, . . . iV — 1 . 

_/V 

M.1: (1 + 2cos0i)(1 + 2cos0 2 )-1 
v.N.2 : (1 + 2cos0i)(l + 2cos0 2 ) - 1 + 2cos20i + 2cos20 2 
M.2 : (1 + 2cos0i + 2cos20i)(l + 2cos0 2 + 2cos20 2 ) - 1 



A discrete Fourier ansatz 

MfcLfcafaij 3 ^) = exp ^^(/ciXi + A; 2 x 2 ) 

yields the same eigenvalues. Setting ki = k 2 = 0, we verify that the maximal 
eigenvalue of A is always equal to the number of neighbors. 

Assuming a large lattice size N, the minimal eigenvalues are determined by 
differentiating with respect to 0,. For the range of eigenvalues A of V~ x A we finally 
obtain: 

v.N. 1: -1 < A fc < 1 
M. 1 : ~ < X k < 1 
13 

v. N. 2 : < A fc < 1 

36 ~ 

29 

M.2: - — <Afc<l. 
96 

(In general, the lower bound given in these expressions refers to the continuous limit 
N — > 00 of and is only approximately realized for discrete 0, i.e. finite N.) For 
the case of global coupling one finds — Af2 1 _ 1 < A < 1 and, in the limit iV — > 00, 
< A < 1. 

The spectrum of a regular graph can be similarly obtained for higher dimensions 
as well as for other lattice types. As an example we mention a triangular lattice (six 
nearest neighbors) which can be considered as a square lattice with the diagonals in 
one direction added. The adjacency matrix is given by 

Arianglc = P + P\ 1 + p2 + P 2 1 + -Pl^2 + Pi ^2 1 J 
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and the corresponding eigenvalues of V 1 A are: 

hiM = 7 ( 2 cos 0i + 2 cos 02 + 2 cos(</>i + <p 2 )) ■ 
o 

The maximal eigenvalue is A max = +1 and the minimal eigenvalue is A m i n = — |. 
This is identical to the case of a square lattice with a 1st order Moore neighborhood. 



3 Stability Analysis for Constant Solutions with 
Synchronous Updating 

Now we consider the stability of a constant solution, 

u(n, x) = u for all n, x , 

of a CML (eq. Q). Given our normalization of parameters, an example of such 
a constant solution is the critical point of the logistic map: u = r — k (As will 
be shown in Sec. 3.3, this is - apart from u — - the only constant solution if 
7 = 0.) Furthermore, we consider a sychronous update of the lattice, which strictly 
corresponds to eq. (JXJ) . Special features related to asynchronous updating will be 
discussed in section HJ 

Let {u k (x)} be a complete set of solutions of the eigenvalue equation for the 
normalized adjacency matrix V~ l A: 

^{V^A) (x, y)u k (y) = — ^ u k (y) = \ k u k {x) . 

y x y~x 

For the stability analysis we consider fluctuations around the constant solution, 

u n [x) = u + 5a n u k (x) , 

and obtain: 

5a k n+1 = (1 - e)f(u)6a k n + e\ k g'(u)Sa k n . (3) 

Hence, the constant solution is stable if e satisfies the following inequalities: 

- 1< (1 - e)f(u) + e\ k g'(u) < +1 for all k . (4) 

Inserting the logistic map for /, and using g as defined in eq. (j2J) we find: 

f\u)=2-r and g\u) = a - (r - 2)0 + (r - 2) 2 7 . 

This leads to the stability conditions: 

r - 3 < e (r - 2 + X k [a - (r - 2)/3+ (r - 2) 2 7 ]) < r - 1 for all k . (5) 

Note that a — (r — 2)f3 + (r — 2) 2 j can be positive or negative. Depending on this sign, 
it is the smallest or the largest eigenvalue of the normalized adjacency matrix for 
which one of the inequality conditions in eq. © is first violated. 

Due to the many parameters (a, (3, 7, e and r) in the CML (eq. Q), and in the 
stability conditions (eq. EJ), we restrict our discussion to special cases emphasizing 
those parameter ranges which have been numerically investigated in previous work. 
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3.1 Pure Causal Coupling 

Pure causal coupling is characterized by (3 = 7 = and a — 1. The stability 
conditions (eq. EJ) then read: 

r-3<e(r-2 + A fc )<r-l. 

The lower inequality defines the critical coupling strength: 

r — 3 



T 2 -f- A m ; n 

Inserting the different minimal eigenvalues for the various neighborhoods we find: 

r — 3 



(6) 



v. N. 1 e c 

M.l e c 

v. N. 2 e c 

M.2 e c 

global e c 



r-2-1 

r — 3 

r - 2.5 

r — 3 r — 3 



13 

36 

r — 3 r — 3 



r-2-±f r- 2.361... 

3d 



_ 29 

r — 3 



r-2-g r- 2.302... 

96 



r - 2 



These results can be compared with numerical simulations presented in |Hj . Fig- 
ure 1 gives a compact representation of those results, showing the critical coupling 
strength e c , beyond which the unstable fixed point of u is stabilized, as a function of 
the control parameter r of the logistic map for different neighborhoods as indicated 
in the figure and for synchronous updating. The case of global coupling turns out 
to be realized for asynchronous updating (cf. Sec. 4) as well. 

Note that for the first order von Neumann neighborhood the constant solution 
is always unstable for e < 1. For r = 4 (right hand side of Fig. 1) we obtain the 
following critical values of e: 



v.N.l 


e c 


= 1 


M.l 


e c 


= 0.66... 


v.N. 2 


e c 


= 0.61... 


M.2 




96 


e c 








163 


global 


e c 


= 0.5 . 



0.589... 



It is easy to check that the numerical results of Fig. 1 agree perfectly well with 
the theoretical results according to eq. (6). 
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0.7 




control parameter r 



Figure 1: Critical coupling strength for stabilization onset, e c , as a function of r, 
3 < r < 4. The individual values are numerical results for synchronous updating. All 
curves are given by the functions of e c on r according to eq. (6) . 



3.2 Causal and Non-Causal Coupling 

Next we consider the case 7 = 0, /? = 1 — a, and r = 4. This situation has been 
investigated numerically in j^j. From eq. (5), we now obtain the stability conditions: 

1 < e(2 + A fe (3a-2)) < 3 . (7) 

For a — |, this becomes independent of the neighborhood. In this case the 
critical value e c derived from the lower bound of eq. (7) is: 

e c (a = §) = -. 

For a > |, the minimal eigenvalue of the adjacency matrix determines the sta- 
bility of the CML with respect to the lower bound in eq. (7) and we find: 

v. N. 1 e ! 



M.l 



4 -3a 
1 

3 - fa 



v. N.2 
M.2 
global 

For a < |, the lower bound is first saturated for the maximal eigenvalue A max = 1 
of the normalized adjacency matrix. Since this is independent of the neighborhood 
and holds for any simple graph, we have the result: 

1 

° 3a 

As e < 1 there is no critical value for e for a < |. In this regime the constant 
solution is generally unstable (for synchronous updating). 

Figure 2 shows the numerical results for combinations of causal and non-causal 
coupling for both synchronous and asynchronous updating. All results for syn- 
chronous updating agree perfectly with our theoretical results. Asynchronous up- 
dating will be discussed separately in Sec. 4. 



36 

6c ~ 98 - 39a 
96 

Cc ~ 250 - 87a 
1 

£c= 2- 



3.3 Causal, Non-Causal, and Anti-Causal Coupling 

This case may be of less practical interest but can shed some light on general features 
of CMLs. In particular, it presents features which are not shared by the coupling 
scenarios discussed in the preceding sections. We distinguish three cases: pure anti- 
causal coupling (7 = 1, a — ft — 0), combinations of causal and anti-causal coupling 
(ft = 0, a = 1 — 7), and combinations of non-causal and anti-causal coupling (a = 0, 
ft = 1 — 7) . Our discussion refers mostly to the stability of the constant solution at 
u = For 7 7^ there are additional constant solutions which depend on e. 

3.3.1 Pure Anti-Causal Coupling: 7 = 1, a = ft = 

The stability conditions in this case are given by: 

r-3<e(r-2 + X k {r - 2) 2 ) < r - 1 . 

For the lower bound, we only obtain solutions in the range < e < 1 if: 

1 



\r-2f 



For r = 4, this condition is not fulfilled for the neighborhoods considered in this 
paper except for global coupling. 

For a given graph, eq. (JBJ) is satisfied if 



r < 2 + 
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Figure 2: Critical coupling strength e c for stabilization onset as a function of the 
degree a of causal coupling for r = 4. Different symbols refer to different neighbor- 
hoods as explained in the figure. The solid line at e c = 0.5 characterizes the situation 
for asynchronous updating, independent of the neighborhood considered. 



In this case, the stability conditions are given by 

r — 3 r — 1 

(r-2)-|A min |(r-2)2 < 6 < (r - 2) + (r - 2f ' 

3.3.2 Causal and Anti-Causal Coupling: f3 — 0, a = 1 — 7 

For r = 4, the stability conditions now assume the form: 

Ke(2 + A fc (l + 3 7 )) <3 . 
There is a general upper bound 

1 



11 



beyond which the constant solution at the unstable fixed point remains unstable. 
The lower bound for e depends on the type of neighborhood: 



v. N. 1 e > — - — nowhere satisfied 

1 — 37 

2 1 

M. 1 e > — if 7 < - 

3-37 3 

36 „ 23 

v. N. 2 e > if 7 < — 

59 - 397 39 

96 67 
M. 2 e > — — — if 7 < — 

163 - 877 87 

global e > - for all 7 . 

Fig. 3 shows numerically obtained CML state values as a function of coupling 
strength e for selected values of 7 and for a first order von Neumann neighborhood. 
There is no stable solution at the unstable fixed point, corresponding to our theo- 
retical result that the stability condition is nowhere satisfied in this case. However, 
we observe e-dependent solutions different from the unstable fixed point. This is 
due to the fact that for constant solutions and 7 7^ eq. (1) is of fourth order. The 
solutions can be characterized (for r = 4) by 



(9) 



Real solutions exist only for e > as can be recognized by the stable regimes in 
Fig. 3b-e. In Fig. 3a there is only a bistable solution. 




3.3.3 Non-Causal and Anti-Causal Coupling: a — 0, (3=1 — 7 

For several interesting values of 7, Fig. 4 shows numerically obtained state values 
of the CML as a function of coupling strength e for a first order von Neumann 
neighborhood, Fig. 5 shows the same for a first order Moore neighborhood. 
The stability conditions for r = 4 read: 

1 < e(2 + A fc (6 7 -2)) < 3 . 

For 7 = |, we find e c = \ from the lower bound, while the upper bound gives 
no restriction. For | < 7 < |, only the maximal eigenvalue is relevant and the 
conditions are independent of the considered graph. We find: 



1 




This is confirmed by the numerical results for a first order von Neumann neighbor- 
hood in Fig. 4b-e and for a first order Moore neighborhood in Fig. 5e. 
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Figure 3: Stability diagram for synchronously updated CMLs with a von Neumann 
neighborhood of order 1. For selected values of 7 — 1 — a, the distribution of states is 
plotted versus coupling strength e. Simulations are based on random initial conditions 
on a 50 x 50 lattice, at least 5000 iterations, and r — 4. 



There is no stability for 7 < |. (The upper bound leads to no restriction.) This 
can be seen in Fig. 4e-f for a first order von Neumann neighborhood and in Fig. 5f 
for a first order Moore neighborhood. 
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Figure 4: Stability diagram for synchronously updated CMLs with a von Neumann 
neighborhood of order 1. For selected values of 7 = 1 — /?, the distribution of states 
is plotted versus coupling strength e. Simulation details are as in Fig. 3. 



For 7 > I we find: 
v. N. 1 

M.l 

v. N.2 

M.2 

global 



1 



< e 



4-67 

1 1 

3-37 < 6 < 27 

18 1 
< e < — 

49 - 397 27 

48 1 

125 - 877 < 6 < 27 

1 1 
— < e < -— 

2 27 

14 
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Figure 5: Stability diagram for synchronously updated CMLs with a Moore neigh- 
borhood of order 1. For selected values of 7 = 1 — f3, the distribution of states is 
plotted versus coupling strength e. Simulation details are as in Fig. 3. 



For the regime | < 7 < |, the theoretical results are confirmed by the example 
in Fig. 4a for a first order von Neumann neighborhood and in Fig. 5d for a first 
order Moore neighborhood. 

While a first order von Neumann neighborhood produces no stable solution for 
7 > 0.5, this is different for a first order Moore neighborhood (and for higher order 
neighborhoods). Figure 5c shows e c = 0.67 at 7 = 0.5 for a first order Moore 
neighborhood. In addition to the theoretically predicted stability of the unstable 
fixed point for all e c < 1, anti-causal coupling (7 7^ 0) yields e-dependent stable 
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solutions as well. For r = 4 they are now given by 

5 1 / 4~ . . 

Real solutions exist only for e > <^ and can be recognized in Fig. 5a-c. In particular, 

Fig. 5c shows a discontinuous transition of the constant solution from u = ^ = f 
to u = | at e = |, which represents a second critical coupling value at 7 = 0.5 in 
Fig. 6. 

The constant solution for a first order Moore neighborhood becomes unstable 
when 7 > I and e < ^. This is confirmed in Fig. 5a-b. Upper and lower bounds 
coincide for 7 = 0.6. 
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anti-causal coupling degree 7 

Figure 6: Critical coupling strength e c for stabilization onset as a function of the 
degree 7 = 1 — (3 of anti-causal versus non-causal coupling for r = 4 and for syn- 
chronous updating. The solid lines represent theoretical predictions and the dots 
refer to numerically obtained values for von Neumann and Moore neighborhods of 
first order. 

All results for combinations of non-causal and anti-causal coupling with syn- 
chronous update for first order neighborhoods are compactly summarized in Figure 
6, representing e c as a function of 7 = 1— (3 (for r = 4). Most numerical data and the- 
oretical curves agree perfectly. Results around 7 = 0.5, deviating from the curves in 
Fig. 6, do still satisfy the stability inequalities. Possible reasons for these deviations 
are due to coexisting stable solutions. This is presently under investigation. 
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4 Asynchronous Updating 



Our results so far were derived for a synchronous updating scenario for the CML 
configurations. For asynchronous updating the analytic treatment is more elaborate. 
The possibility of different types of asynchronous updates makes this even worse. 

For synchronous updating, we can interprete the CML (eq. as the iteration of 
a single mapping F which maps a vector {u(n, x)} (with components x) to a vector 
{u(n + 1, x)}: 

F : {u(n,x)} — > {u(n+l,x)} 

where u(n + l,x) is given by eq. (0) for all x (i.e., synchronously). The linear 
approximation around the constant solution (leading to eq. Oj) can be written as 

5u n +i(x) = y j F'(x,y)8u n (y) 
v 

where the matrix F'(x,y) is given by: 

F'(x, y) = (1 - e)f'(u)8(x, y) + eg'(u)(V- l A)(x, y) 

with 5(x,y) as Kronecker's 5. The stability condition (eq. 0} is simply a condition 
for the eigenvalues of F' . 

By contrast, asynchronous updating corresponds to a set of mappings (one for 
each vertex z) such that: 

F z :{u(n,x)} — ► {u(n + l,x)} , 

where now {u(n + 1, x)} is given by: 

u(n, x) for x ^ z 
"(" + !■•'■) = < (l-e)/( u (n,x)) + — 5^<7(ti(n,x)) for x = z ■ 



Correspondingly, the linear approximation now depends on the vertex z, for which 
the update has been made: 



F' z {x,y) = { 3 X,V \ T X ^ Z 

z\ mi | F'(x,y) for x = z 



The stability condition is now a condition for the eigenvalues of a matrix, which is a 
product of F^-matrices, where the product runs over z either randomly or according 
to some specific sequential rule. In general, the matrices F' z (x,y) do not commute 
for different values of z which makes it difficult to determine the eigenvalues of this 
product. 

However, despite this difficulty to derive an exact analytical solution for the 
stability of asynchronously updated configurations, there is an intuitive argument 
for the stability of such configurations. 
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Consider some perturbation Su(x) around a constant solution u of the CML. 
The lower bound of the stability equation (eq. Ej) corresponds to a change of Su(x), 
including a change in sign, at the vertex x where the update is made. In particular, a 
perturbation proportional to the eigenfunction of the adjacency matrix for which the 
corresponding eigenvalue violates the lower stability bound for the constant solution 
is simply multiplied by —1 at the critical value of e. 

When an update is made at some vertex z, we can assume that its neighbors 
have been randomly updated before and therefore their signs have changed randomly 
This implies that, on average, the neighbors of vertex z cancel and do not contribute 
to the update, so we have an effective update equation of the form: 

Su(n + 1, z) = (1 — e)f'(u)8u(n, z) . 

The lower bound of the stability condition of this equation reads: 

-1< (1 - e)f{u) , 

which for the logistic mapping becomes: 

r — 3 



This equation has been derived from different arguments in [§]. It corresponds to 
the solid curve shown in Fig. 1. 

However, if e approaches the upper bound of the stability inequality, a per- 
turbation 5u(x) does not change its sign. Therefore, even if a random number of 
neighbors of a vertex z, which is to be updated, have already been updated before, 
their fluctuations will not necessarily cancel on average. This situation is similar to 
synchronous updating. This argument makes it plausible that the upper bound for 
e, as derived from synchronous updating, is also valid for asynchronous updating. 

Among the various situations studied numerically in 0111], there is only one case 
where the upper bound of the stability inequality becomes relevant: combinations 
of causal and non-causal coupling in the regime a < | for a first order von Neumann 
neighborhood. For this case the upper bound stability condition ((7j) reads: 

3 



This has been numerically confirmed for a number of examples shown in Fig. 2. 

5 Generalized Neighborhoods 

In the preceding sections, all neighbors of a vertex have been treated identically. 
However, it is possible to consider more general, inhomogeneous couplings of the 
vertices of the graph on which the CML is defined. For instance, the coupling 
strength may depend on the distance on the graph (where distance referes to the 
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length of the minimal path). Such generalized couplings are also suggested by renor- 
malization methods, which, even if the starting point is a nearest neighbor coupling, 
effectively lead to couplings between more distant points. 

The most general case is given by a symmetric adjacency matrix with a non- 
negative coupling between any two vertices: 

A {9) (x,y) = A i9) (y,x)>0 A ia) (x, x) = . 

We can define the generalized valence matrix by: 

V (g) (r y) = / £^ (ff) (^) for x = y 
y ,y) \ otherwise ' 

and we assume that V^ 9 \x, x) > for all x. As in Sec. 2, one can show that 
(1 ± A^) is non-negative and that the largest eigenvalue of A^> is +1 

(corresponding to the constant function on the graph). In the context of Markov 
processes the matrix A^ is called a Markov matrix. 

A simple example is given by a mixture of von Neumann and Moore neighbor- 
hoods of first order, where nearest horizontal and vertical neighbors are linked by a 
coupling ai and the additional diagonal neighbors are linked by a coupling a 2 . We 
now define the generalized adjacency matrix 

{a\ if x and y are nearest neighbors 
a 2 if x and y are diagonal neighbors . 
otherwise 

This matrix can be written in the form: 

A(x, y) = ai(Pi + Pf 1 + P 2 + P 2 - 2 ) + a 2 (pP 2 + P X P 2 ~ X + P{ l P 2 + P^P 2 l ) , 

and the corresponding eigenvalues A of A are: 

A = 2ai(cos 0i + cos 2 ) + 4a 2 cos 0i cos 2 . 

The maximal eigenvalue (equal to the generalized valence of a vertex) and the min- 
imal eigenvalue of A are: 

Amax = v = 4(ai + a 2 ) and A min = -4max(ai, a 2 ) . 

The normalized minimal eigenvalue can now be tuned between the case of first order 
von Neumann (A min = — 1) for a 2 = (equivalently, a x = 0) and first order Moore 
(A min = ~\) for ai = a 2 . 

An interesting case occurs if the coupling a(L) depends only on the lattice- 
distance L between two points such that a(Li + L 2 ) = a(Li)a(L 2 ). In this case 
a(L) = p L for some positive constant p which we assume to be smaller than 1 in 
order to avoid divergencies on infinite lattices. For the following we assume the 
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lattice to be rectangular, <i-dimensional and consider the limit of infinite lattice size, 
N — > oo. The adjacency matrix can then be written in the form: 

d / oo \ 

a = U[ 1 + J2p l ^ l + p ^)- 1 

i=l V L=l / 

with eigenvalues: 

d / oo \ d , 1 2 \ 

For the maximal and minimal normalized eigenvalues we finally obtain: 

A -1 A - ( l - p 

\1+P 

The significance of this exercise is obvious. Already for the 1-dimensional lattice 
(the circular graph) we can mimick, by a suitable tuning of p, the largest and smallest 
eigenvalues of any other arbitrary graph or generalized coupling. These eigenvalues 
are independent of the functions f(x) (for which we used the logistic map) and 
g(x). The stability of the constant solution of a CML does not depend on any other 
property of the graph except these two eigenvalues. 




6 Summary 

Coupled map lattices (CMLs) and other complex model systems produce an over- 
whelming phenomenological variety of features which are notoriously difficult to 
classify and explain in a compact and comprehensive way. A promising concept for 
such purposes is the concept of stability since it characterizes complex systems in 
their different dynamical regimes in a most basic fashion. 

The stability of CMLs has been studied numerically by a number of authors 
before. A large class of such studies were carried out in the context of synchro- 
nizing and controlling complex systems. In contrast, analytical tools explaining the 
numerical results have not been broadly applied so far. Insofar as CMLs can be con- 
sidered as implemented on special types of graphs, spectral graph theory provides a 
particularly convenient option to analyze the stability of CMLs theoretically. The 
maximal and minimal eigenvalues of the normalized adjacency matrix of the graph, 
on which the CML is defined, determine its stability properties. 

We used these tools to analyze the stability properties of CMLs for different types 
of neighborhoods, different coupling scenarios and for different updating rules. In 
this paper we focused on the stability at or around temporally constant solutions of 
the CML as a whole. Particularly interesting is the case of globally stable solutions 
at locally unstable fixed points, which are independent of the coupling strength e. 
In addition we discovered constant solutions which are explicitly e-dependent. 
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For synchronous updating of the CML, we derived exact stability conditions for 
causal, non-causal, and anti-causal coupling. These different coupling scenarios dis- 
tinguish whether states of a vertex of the graph are assumed to interact with their 
preceding, simultaneous, or future neighboring states, respectively. The considered 
neighborhoods are of low-order von Neumann and Moore type. All theoretical re- 
sults agree excellently with numerical studies published earlier. Novel numerical 
simulations, referring to situations not treated before, additionally confirm our the- 
oretical results. 

For asynchronous updating, the stability analysis is in general more involved than 
for synchronous updating. In the case of a random updating sequence, the derivation 
simplifies and, again, gives results that agree with those of earlier numerical studies. 
A heuristic stability condition derived earlier could be confirmed as well. 

Finally, we generalized our approach to inhomogeneous neighborhoods, such as 
mixtures of neighborhoods of different types or coupling strenghts depending on the 
distance from a vertex. It could be shown that the stability of constant solutions 
of CMLs is completely characterized by the maximal and minimal eigenvalues of 
its normalized adjacency matrix, independent of the dynamics of the vertices and 
independent of the causality features assumed for coupling. 
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